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Abstract 

In this paper a certain class of singular integral equations that may 
arise from the mixed boundary value problems in nonhomogeneous materials 
is considered. The distinguishing feature of these equations is that in 
addition to the Cauchy singularity, the kernels contain terms that are sin- 
gular only at the end points. In the form of the singular integral equations 
adopted, the density function is a potential or a displacement_and conse- 
quently the kernel has strong singularities of the form (t-x)~S x n “ 2 (t+x) n , 
(n^2, 0<x , t<b ). The complex function theory is used to determine the funda- 
mental function of the problem for the general case and a simple numerical 
technique is described to solve the integral equation. Two examples from 
the theory of elasticity are then considered to show the application of the 
technique. 


1. Introduction 

In elasticity or in potential theory if the medium contains a planar 
imperfection representing a discontinuity in displacements or potential (e.g., 
a crack, a plane insulation, a barrier), or in stress or flux vector (e.g., 
a plane inclusion, a distributed source), the related mixed boundary value 
problem may be formulated in terms of either a system of dual integral equa- 
tions [1] or a singular integral equation [2,3]. If the "cut" corresponding 
to the plane of potential or flux discontinuity is fully embedded in a homo- 
geneous component in the medium, it is known that the dominant part of the 
integral equation contains only a simple Cauchy kernel associated with the 
corresponding infinite domain and the remaining geometry of the medium is 
represented by a Fredholm kernel . On the other hand, if the medium is 
nonhomogeneous with discontinuous material parameters and if the cut inter- 
sects such a plane of discontinuity, then it is also known that the dominant 
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kernel of the integral equation contains other singular terms in addition to the 
Cauchy kernel [3], These kernels (which are sometimes known as the general- 
ized Cauchy kernels) become unbounded as the variables approach the end 
of the cut responding to the point of intersection with the bimaterial inter- 
face. 

In the case of singular integral equations with ordinary Cauchy kernel 
the solution of the problem is quite straightforward and there are very highly 
efficient numerical techniques to obtain it. Most of these simple techniques 
are based on a variety of quadrature formulas developed for singular inte- 
grals (see, for example [4]). Another such technique introducing further 
computational economy as well as maintaining the basic numerical simplicity 
in the applications was discussed in a recent article [5], The main features 
of this technique are (i) potential rather than flux type quantities are 
used as the unknown functions in deriving the integral equations resulting 
in a dominant kernel with a singularity stronger than that of the standard 
Cauchy kernel, (ii) the unknown function is expressed in terms of the fun- 
damental solution of the problem and a bounded function containing a set 
of unknown coefficients, (iii) interpreting the integrals with strong singu- 
larities in the Hadamard sense, some useful formulas are developed to 
evaluate the integrals corresponding to the dominant part of the integral 
equation in closed form, and (iv) an appropriate collocation method is used 
to solve the resulting functional equation for the unknown coefficients. 

In [5] this method was used to solve also a special case of singular integral 
equations with a generalized Cauchy kernel, namely that of a cut intersect- 
ing a free boundary, and was shown to have distinct advantages over the 
conventional quadrature methods with regard to accuracy and computer time. 

In this paper the mixed boundary value problems leading to singular 
integral equations with a generalized Cauchy kernel are reconsidered by 
formulating them in terms of "potentials" as the unknown functions. The 
terms in the resulting dominant kernel would then have strong singularities 
and must be interpreted in the Hadamard sense. It is shown that the complex 
function theory can again be used to determine the fundamental solution. 

The main objective of the paper is to develop a simple and efficient tech- 
nique for solving singular integral equations with generalized Cauchy kernels. 
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2. Integral Equations 

Consider first the simple problem in potential theory described in Fig. 
1. In the terminology of elasticity the problem is one of antiplane shear 
loading of a nonhomogeneous medium which consists of two bonded half spaces 



Fig. 1 Bonded Half-Planes (Anti plane Shear Loading) 

with shear moduli u-j and y 2 * Tire medium contains a crack along 0 =tt, a<r<b 
and, in the perturbation problem considered, the crack surface traction 
CTle z (r,i r )=q(r) is the only nonzero external load. If w-j and w 2 are the z- 
component of the displacement vector in materials 1 and 2, respectively, it 
was shown that (see, for example, [6,7]) the differential equations and 
boundary conditions 


V 2 w-J = 0, v 2 w 2 = 0 , 

(la,b) 

W 1 = w 2* ° 19 Z = °2ez * (0 = i ’ 0 - r< “ ) 

(2a,b) 

w 2 (r,0) = 0,. (0<r<») 

(3) 

w-j(r,ir) = 0, (0<r<a, b<r<») , 

(4a) 

cr i qz ( r ,7r ) = q( r ) > (a<r<b) 

(4b) 
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would give the following integral equation: 


b 

J-J- (Fr + Fr^ ^( t )dt “ ~ q(r) , (a<r<b) (5) 

a 

where 

^(r) = [w-j(r s ir+0) - w-j(r,7T-0)] , (6) 

X = (^i )/ Cvt-j + ii2 ^ (7) 


and, because of symmetry, only one half of the medium is considered. We 
note that for a=0 and y-|^p 2 the term *(t+r)“^ becomes unbounded for r=0=t 
and the kernel is one of generalized Cauchy type. 

Integrating (5) by parts and using (4a) it may be seen that 
b 

7 f + it7FF Mt ) dt = ^ q(r) • (a<r<b) • (8) 

a 

where 


W(t) = W-|(t,TT+0) - W-j ( t , TT**0 ) . 


(9) 


The integrals in (5) and (8) are to be interpreted in Cauchy principal value 
and Hadamard sense, respectively [5], The integral equation can also be 
obtained by using w(t) rather than ^(t) as the unknown function and, for 
example, following the procedure outlined in £6], It may be noted that in 
terms of w some of the stress components of physical interest may be expressed 
as 


a 2ez 


(r,0) = - 


Vt-iU 


1 2 
M,+« 2 


1 

TT 


w(t) 

Tt+F) 


T dt , (0<r<») 


(10) 


a lez (r » w/2) = a 2 ez {r>7r/2) 


¥2_1 ' 
Pl+P 2 TT . 


w(t)dt, (0<r<») , (11) 


b 

ff lez( r,ir ) = 2 7 | C (t- r) 2 * ~ (tFr ) ( t ) dt » (°l r<a > b<r<») . (12) 

a 
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Similarly, for the crack geometry shown in Fig. 1 under the following 
symmetric in-plane loading condition 


a la0 (r,Tr) = -p(r) , a lrQ (r, ir) = 0, (a<r<b) , 


(13) 


the integral equation may be obtained as (see, for example, [8] for the pro- 
cedure) 

b b 


j dt + | K(r,t)v(t)dt = -ir p(r), (a<r<b) 


where 


v(t) = u l0 (t,Tr+O) - u l0 (t,7T-O), (a<t<b) , 


2c 2 r 3c 3 r2 


K ( r,t ) = Tt+rF + Tt+rF + Tt+rP * 


= I [1 (1+lc l )m _ 3(l-m) 


'1 2 


m+tc 


2 l+mic-j 


] , 


- = 6 l.ljOl c = . 4(1 -ml = / ■ 

'2 1+m^ » c 3 1+mK-j » m V W 1 * 


(14) 

(15) 

(16) 


(17a-d) 


In the above formulation y-j and u 2 are the shear moduli, k.. s 3-4v.. for plane 
strain, K.j=(3-v.j )/ ( 1 +v^ ) for plane stress, (i=l,2), and v-j and v 2 are the 
Poisson's ratio of the two materials. Once the integral equation is solved, 
in this problem too all the desired field quantities may be expressed in 
terms of v(t) and the corresponding kernels. For example, the cleavage 
stresses along the planes e=0, e=rr/2 and e=ir may be written as [8] 

b 

V 2 f rr 3 1 \ 1 , { 1 4r l v(f)dt 

L ^m +<2 1+mic^ Tt+rp" M+m<^ m+K^ (t+r) d ^ ' * 


a 200 (r,O) - ^ 


(0<r<«) 


(18) 


CT 1 0Q ( r » 7T / 2 ) = cr 2 (r,7r/2) = — 


f=r-C 1 


m+i <2 L t z +r 2 


2r 2 -I 
(t z +r 2 ) zJ 


- T^- Eplpr ' T tW + TtW ]>v(t)dt ’ (0<r< “> ’ (19) 
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(r, tt) = 


( 20 ) 


CT 1 00 


2y 1 

TT { 1 +K-j ) 


b 

» 

a 


[f t-r ) 2 + K ( r » t ) 3 v (t)dt, 


( 0 <r<a, b <r«») 


9 


where K(r,t) is given by (16), 

In the problems given by the integral equations ( 8 ) and (14) for a crack 
embedded in a homogeneous phase, (i.e., for a> 0 ), the dominant kernel con- 
sists of (t-r) , the fundamental function of the integral equations is 

(b-r) 2 (r-a)% and the solutions may be obtained by using the technique 
described in [5], On the other hand, if a=0 clearly the kernels are of the 
generalized Cauchy type and before attempting to solve the problem the fun- 
damental function of the integral equation needs to be determined. This may 
easily be done by using certain properties Cauchy integrals [2]. 


3. The Fundamental Function 

For a=0 we note that the general form of the integral equations ( 8 ) 
and (14) is 

r K B x k ‘ 2 

“f [t v y -r ? + £ k “ 2 k ] f(t)dt = p(x), ( 0 <x<b) . ( 21 ) 

T (t-x) k =2 (t+x) k 

To determine the behavior of f(t) near and at the end points, following [2] 
we def i ne 

f(t) = g(t)t 6 (b-t) a = g(t)<r”V(t-b) a . 

(Re(a,e)>0, g(0)ji0, g(b)!f0, 0<t<b), (22) 

B] cc| 

where g(t) is bounded in the closed interval [ 0 ,b] and the function t (b-t) 
constituting the leading term is generally known as the fundamental function 

of the integral equation (a^^,..., B=B-|,B 2 >»**» 0 <Re(ct^ )<Re(a 2 ) 

0<Re ( B-j ) <Re ( 3 2 )>•••) • 

Defining the sectional ly holomorphic function 
b 

F(z) = 7 J dt » (z=x+iy) » (23) 

o 

it may be seen that 
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(24) 


1 1 dt = G(z > = 3? F « = ^ \ dt ' 


From (24) the asymptotic behavior of G(z) near the ends may be expressed as 

G(z) - -9(0)b« egg. z*-' + 9 <b)b S ^ (z-bf* 1 ♦ G 0 (z) , (25) 

where G Q (z) is bounded everywhere except possibly the end points near which 
it may have the following behavior 

P •-! 

l G 0 ( z )| <c il z ~ d ii 1 » ( d 1 =0> p-,>Re(s)v d 2 =b, p 2 >Re(a)). (26) 

First we note that if z is not on the cut, then G(z) is holomorphic and 
at a point z=z o =xe 10 , (0<x<b), we can write 


7 1 T dt = G(z o> = - 9(0)b “ ~ 


i[(e-l) 0 -irB] 
sirnrg " 


where G^ is bounded- everywhere and near z Q =0 has a behavior similar to that 
of G Q . For example, if 9=ir (i.e., if z Q =-x), by differentiating from (27) 
it may easily be shown that near the end x=0 we have 


7 I 


p,+l-k 

+ 0(x ) , (k>2, x>0, p 1 >Re( e) ) . 


Similarly, substituting e=ir/2 and e=-ir/2, (i.e., for z Q =+ix) and using 
i = V ( 1 _ l \ t 1 / 1 . 1 \ 

t z +x 2 2ix l t-ix " t+TX J ’ t 2 +x z ~ 1 't-ix t+i x } 


it can be shown that 
b 


U ? _1 
1 f dt =— Tb g(0)b a x s " 1 +0(x 1 ), (x>0, p-| >Re{ s) ) , 


we find 


b 

' { dt * l&llL t_3)...( R +,.,, l 

(MjitoTf m 9(0)bV+l-3k 

. Pj+!-2fc 

+ 0 (x 1 ) f. „ 

» i 2 » *>0, p^Refs}} 

+ n/ P ? +2 ‘ 2k 

J ' ( k i2, x>0, Pj>Re(g)) . 

% u S f„ 9 the Plemelj f„ mula [2] 

'( (0<x<b) 

nd 


,« x e+2-2k 


* / TE^x|r dt = J. ^ f f( t ) 

0 ^V^ dt = r^W Wj>(0<x<b) 

";’:r ,e>us,n9p --— 

t-2 7 = r. _ - 


' ' V t CJi Q7- ■ 

r,8-v , ' 

1 J a fx e 70 J^I . «-] 
e->+0 x , 

s Tim r x gie-, 

9-^-0 J 


'V' 1 = X s-I e «2* e 


fr ™ (34) and (25) 

b ' n f °Hows that 

I 1 f[t) . 

* J o n^r * * -9(0Jb« s cotM x e-’ 

' Sfb > be “C0t(™ >f5 . x) c-l + 0( /T-' >(b _ x) P 2 - I) 

( °<x<b, P,>R e (a), p2>ReW) 


(31) 


(32; 


(33 


(34 
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Substituting now from (36) and (28) into the integral equation (21) we 
obtain 


-g(0)b c ‘x B " : '{8cot(Tre) 


K 

E 

k=2 


k-2 


(- 1 )‘ 


B(8-1) v .(e-k+2) 1 
(k-1 ) l s i n 6 i 


Pi -1 Po-1 

-g(b)b a acot(Tra) (b-x) a +0(x , (b-x) L ) = p(x), (0<x<b) . 


(37) 


Now, by observing that g(0)j f 0, g(b)^0, p 1 >Re(B)>0, p 2 >Re(a)>0 and p(x) is 
bounded in 0<x<b, from (37) multiplying both sides first by x 1-B and letting 

1 -a 

x->0 then by (b-x) and letting x-*b, for the leading terms 0<Re(ct,8)<l we 
obtain the following characteristic equations: 


K \k 

BCOtirB - E B. ' 

k=2 


B(8-l)...(B-k+2) _ n 
"(k-t) isiniTg " 0 


COtira = 0 


(38a, b) 


Note that for the end x=b that is embedded in a homogeneous component of the 
medium (38b) gives the known result a =1/2, 

For the integral equation (8), if a=0 by defining 

w(t) = g(t)t 0 (b-t) a , (0<t<b, Re(ot,8)>0) , (39) 

from (36) and (37) we find 

COS it 8 - X = 0 , cot TO = 0 . (40a ,b) 

Similarly, for a=0 in the in-plane elasticity problem given by (14) if we 
let 

v(t) = g(t)t B (b-t) a , (0<t<b, Re(a,8)>0) (41 ) 

from (14), (16), (21) and (38a, b) it may be shown that 
-cosire + c^ + c 2 (l-e) + c 3 (1-8)(2-b)/2 = 0 , 

(42a, b) 

COtira = 0 
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We also note that after determining the unknown function g(t), the aysmptotic 
behavior of the stress components around the end points x=0 and x=b may be 
obtained by substituting from the asymptotic relations such as (28) -(30) 
and (36) into the expressions of stresses (e.g., (9)-(ll) and (18)-(20)). 

4. The Numerical Technique 

To solve the integral equation (21) having the generalized Cauchy kernel 
we express the unknown function as 

f(t) - g(t)t S Vt) % , (43) 

where ^ is the smallest positive root of (38a)and g(t) is an unknown bounded 
function. By defining the following normalized quantities 

(44) 


(45) 

*(t) = h(T)(l-T)*(l+T) , (-1<T<1) , (46) 

where h(x) is the new unknown function. As in [5], the numerical solution 
of (45) may now be obtained by expanding h(r) into a series with known coor- 
dinate functions and unknown coefficients and by reducing the resulting func- 
tional equation to a system of algebraic equations through a suitable collo- 
cation technique. No particular system of coordinate functions seems to 
have a special advantage and a simple power series such as 

N 

h(x) = i a x , (— 1 <t<1 ) , (47) 

n=0 n 

seems to be quite adequate. Also, expressing the finite-part integral as 
[5] 


t = | (1+t), x = | (1+p), f(t) = | 4 >(t), p(x) = p-, (p ) , 
from (21) we obtain 

if 1 ^ ^k-2^** 3 ^ 

— T [ / " p \2 * Z — 1 V — — ]<f>(x)dT = P-j(p), (-1 <p<l ) , 

" i. (T " pJ k=2 (r+ P +2) k 1 
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1 


-1 




[4'(t)-4>(p)-(x-p)4>*(p)] 


-1 


dr 

( T ~p) 


- *(p)(t^ + tV } + 


(48) 


(*) 


seems to give very satisfactory results^ ' . Higher accuracy is obtained in 

the numerical solution if the density of collocation points is increased near 

the ends by, for example, selecting p . ,(j=0,l ,N) as the roots of Chebychev polynomials , 

J 

Numerically the technique described above would give fast converging 
results provided 


a.j ~ot-j 1 , ^ 1 , (i,j 2,3,...) . 


(49) 


On the other hand if the characteristic equations contain roots for which 
a^-a-j <1 or Sj-B-|<l, (i»j=2,3...), then the numerical technique would converge 
slowly and would have to be modified. To see this we note that for the prob- 
lem under consideration formulated in terms a potenti al -type quantity f(t) 
(see (21)) the physically acceptable roots of the characteristic equations 
are ct-j, a 2> ... and » B 2 » - - - wlth Re («-j »6j)>0, (i ,j=l ,2,. . .) and the solu- 
tion may generally be expanded as 


<30 Qjj Ct • 3 * 

f(t) = z z A. .(b-t) \ 3 , (0<t<b) 
i=l j=l 1J 


a U e l 


(50) 


If we now identify the fundamental function as (b-t) t and express 

f(t) = g(t)(b-t) ai t 61 , (51) 


it is seen that 


CO 00 


a r a i + B r e i 


g(t) = z z A.. (b-t) t 
i=l j=l 1J 

Since a.j>a-| and , (i,j=2,3,...) g(t) would be bounded at the ends t=0 


(52) 


FT 


Note that for x->p the integrand on the right hand side of (48) is 
<J>"(p)/2. 
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andt=b. On the other hand if we consider the derivative of g(t) 

g'(t) = LI A- -[(a .-a-. ) (b-t) 1 t 3 '+(6,-0! )( b-t) 1 4 3 1 ], 

i=lj=V 3 11 31 


( 53 ) 


it is seen that if there are characteristic roots for which 


Bj-B-|“l<0 , (i ,j=2,3,.. . ) (54) 

then at the corresponding end point g 1 ( t ) would become unbounded and conse- 
quently g(t) would be ill-defined. Since g(0) and g(b) invariably represent 
the magnitude of parameters of primary physical importance, it would then 
be necessary to modify the numerical procedure in order to compute these quan- 
tities more accurately. 

As will be shown by the examples given in the next section, the conver- 
gence of the calculations can be improved quite considerably by introducing 
the following simple modification. Let (54) be valid for a i , (i=2,3,. . . ,L) 
and Bj, (j s 2,3,. . . ,M). The corresponding terms are then embedded into the 
definition of f(x) as 

a-. B-i L a- Bi M a, ft.* 

f(t) = g-j (t) (b-t) t + E^g il (t)(b-t) t +^g lj (t)(b-t) t 3 . (55) 

The functions g.j(t) are determined by again expanding them, for example, 

* vJ 

into power series with unknown coefficients and by using a suitable colloca- 
tion method. 


5. Examples 

First we consider the mixed boundary value problem in potential theory 
for the nonhomogeneous medium described in Fig. 1 and formulated by the 
integral equation (8). For a>0 the fundamental function of the integral 
equation is (b-t) z (t-a) % the unknown function can be expressed as 

w(t) = g(t)(b-t)^(t-a) % , (a<t<b) , (56) 

and the stress intensity factors or the parameters giving the strength of 
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of the flux singularity at the end points r=a and r=b may be defined and 
evaluated from 


k 3 (a) 


lim /2(a-r) a l0z (r,Tr), (r<a) 


-J- lim -^ 1 — = 4- g(a)/(b-a)/2 , 
2 t-a /ZXwT 2 


( 57 ) 


k 3 (b) = lim /2(r-b) ^ (r,7r) , (r>b) 

r-4) 

= 4- lim g(b)v/(b-a)/2 . 

£ t-»b /2(ft-t) * 

The problem is solved by introducing the normalized quantities 


(58) 


t = ^4 t + , r = P + > w (t) = 4> ( x ) > 4>(T)=h(x)/l - t 2 , (59) 


b+a 


b-a 


b+a 


b-a 


and by letting 
N 


h(x) = E^x" , Pj = cos(^iLf-) , ( j=0,l . ,N) 


(60) 


The integrals that arise in this solution are evaluated by using the 
following formulas [5]: 

.1 

+ _ k 

(x-p) z 


■ 1 


n 

E 

k=0 


rr? dx - e c^p > (n>0) , 


c k = < 


|0 , for n-k = odd 


^+1 1 ' ? ' 

— - - n - y +9 » for n-k = even. 

2/iT r(2^) 


(61a, b) 


For a material pair corresponding to aluminum and epoxy and for a uni- 
form shear stress q(r)=q Q applied to the crack surface some results are 
given in Table 1, where the previous results calculated by using a 
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Gauss-Chebyshev quadrature technique [6] are also displayed. From the table 
it may be observed that the method seems to give rather good results even 
with a relatively small number of terms used for approximating the unknown 


Table 1. Stress intensity factors for a crack near an 
interface (anti-plane shear). 


p 2 

b+a 

N+l 

k 3 (a) 

k 3 (b> 

y l 

b-a 

%]¥ 

n ft"* 

23.077 

2 

4 

0.958968 

0.976071 



Ref. [6] 

0.95897 

0.97607 



6 

0.712470 

0.924184 


1.1 

10 

0.712075 

0.924190 



Ref. [6] 

0.71208 

0.92419 

0.04333 

1.1 

6 

1 .440085 

1.109604 



Ref. [6] 

1.44009 

1.10958 


function h, the convergence improves as the crack moves away from the inter- 
face, and at the crack tip near the interface the convergence for y 2 >y-j 
appears to be slower than the case for which y-]>u 2 * Analytically, the dif- 
ference in these two cases is only in the sign of X in (8) and, for a>0, is 
difficult to relate it to the convergence. However, as mentioned in the pre- 
vious section, for a=0 there is a good reason for slew convergence and phy- 
sically the two processes are clearly related. 

Consider now the case of a=0 for which the characteristic equations 
are given by (40). For various material combinations the first two roots 
of (40a) are shown in Table 2. In this problem 

w(t) = gftKb-tjV 1 , (62) 
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and the stress intensity factors are defined hy 


k 3 (b) = lim t/F(r-b) <J lez (r»iO > (63) 

r-^b 

1 — 3-i 

k 3 (0) = lim y/2 r a 2ez^ r, °^ * (64) 


Table 2. The roots of the characteristic equation 
COSttB-A=0 


m 2^1 

A 

*1 

6 ? =2-b 1 

0.04333 

0.91694 

0.13065 

1.86935 

0.1 

0.81818 

0.19498 

1.80502 

0.25 

0.60000 

0.29517 

1 . 70483 

0.5 

0.33333 

0.39183 

1.60817 

1. 

0.00000 

0.50000 

1.50000 

2. 

-0.33333 

0.60817 

1.39183 

4. 

-0.60000 

0.70483 

1.29517 

10. 

-0.81818 

0.80502 

1.19498 

23,077 

-0.91694 

. 0.86935 

1.13065 

40. 

-0.95122 

0.90017 

1.09983 


Substituting now from (12) and (36) into (63) and from (10) 
(64) and observing that siniT^=24i|^/(v-|’Ht2) we 


y l b 


5 1 


k 3 (b) =^.H_ g ( b ) 


9 


and (28) into 


(65) 


k 3 (0) — Ai-jU 2 e i 9 (0)i/B7T 


( 66 ) 


For q(r) = q Q a closed form solution of the integral equation (5) is 
given in [7] and [9] which, in terms of the normalized quantities 


5 = t/b, ^(0 = 4>(t) , (0<t<b, 0<C<1) 


(67) 


may be expressed as [10] 
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*T ( «)=--7^ ct 


— + 1) + ( — 5 — -)■ '( 


,-Y, 


y^sin l+ZR 7 " /Np - 1+/R 7 /TT Z 


l)], Y=i-e r 
(68) 


By observing that 4>(t)=dw/dt, from (68) the exact expressions of the stress 
intensity factors may be obtained as 

1-6, 


k 3e (b) = - ^rlijn *(t) = -q Q ^-- 3 ( ^ /2) - , 


1 - 6 , 


k 3e (0) = '^-\v 2 n Jim t 1 >(t) = ^-^72 7 • 


6 ] - B l 

v 1 


(69) 


(70) 


The numerical results given in Tables 3-5 have been normalized as follows: 


Mo) 


UO) = T_ o' » k o(h) - 

3 q 0 (b/2)'- e l 3 


M b > 

q 0 ^7^ 


(71a, b) 


Tables 3 and 4 show the calculated results obtained by assuming the 
solution of (8) in the form 

is 6] b N n I, 6, 

w(t) = g(t)(b-t)^t ' = 7 E a n T n (l-T)^(l+T) 1 

n=0 n 


(72) 


and by following the procedure described in Section 4 above with the collo- 
cation points given in (60). The tables also show the exact values of the 
stress intensity factors calculated from (69) and (70). It is seen that 
for y 2 <ui excellent agreement is obtained by using only six terms in (72). 
However, for y 2 >u-j table 2 shows that e 2 - B-j -1 <0 and as indicated in the 
previous section one would expect the convergence to be slow. This may 
indeed be seen from Tables 3 and 4. 
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Table 3. Normalized stress intensity factors for a crack terminating 
at the interface (anti-plane shear) (N+l=6 terms are used 
in the expansion (72)). 


y 2 

k 3 (0) 

C 3< b ) 

W 1 

(Exact) 

(Exact) 

0.04333 

mm 


1.2558 

1 .2558 

0.1 



1.1940 

1.1940 

0.25 

— 

0.3100 

1.1144 

1.1144 

0.5 

— 

0.5575 

1 .0958 

1.0958 

1. 

1.0000 

1 .0000 

1.0000 

1.0000 

2. 

1 .8036 

1.8134 

0.9598 

0.9598 

4. 

3.2694 

3.3558 

0.9334 

0.9334 

10. 

7.0560 

7.8231 

0.9146 

0.9146 

23.077 

13.6113 

17.3673 

0.9067 

0.9067 

40. 

20.3490 

29.6029 

0.9040 

0.9040 


y 2 

Table 4. Stress intensity factor at the interface, — = 23.077. 

W 1 


N+l 

y°> 

6 

13.61 

10 

14.49 

12 

14.76 

16 

15.14 

20 

15.36 

25 

15.59 

Ref. [6] 

13.13 

Ref. [11] 

14.0 


Exact 17.37 
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Referring to Section 4, for U 2 >u-| we now define the unknown function 
w in (8) as follows: 


w 


v i_Dr/T \^2 / 1 \ - _ n./-, \ ■'i t , . \ c. 


N 2 


(t)=g(t)(b-t)^t 1 =| [(1-t)^(1+x) 1 e a n T n +(l-T)^(l+T) 1 z bx n ], 

c n=0 n n =0 n 


(0<t<b, “1 <t<1 ) , 


( 73 ) 


The collocation points are again selected as the roots of Chebyshev polyno- 
mials, namely 


p j = cos( r+fiW f> • <3 =0 * 1 v->< n i +n 2 +1 )) • (74) 

For two values of y 2 /y-| the results obtained by using the expansion given in 
(73) are shown in Table 5. It is seen that without the second terms in (73) 
the convergence is very slow, whereas with the second term included very accur- 
ate results can be obtained even with a few terms in each series. 

As a second example we consider the plane elasticity problem described 
by Fig. 1 and equations (12)-(17). In this case the characteristic equations 
are given by (42a, b) . Examination of (42b) and the results given in 

[8] show that in all material combinations considered b-j is real and * 

No convergence difficulty is therefore expected. The stress intensity fac- 
tors at the crack tips are defined by 

1-Pi 

k-j (0) = lim /I r » (75) 


k.(b) = lim /2(r-b) ct, (r,ir) 
r-*b 


(76) 


By using (18), (20), (41), (28) and (36), from (75) and (76) we obtain 

0n 


ki (°) ~ y 2 s i n tt 3-j ^m+K 2 " T+mK^ 


+ 2(1 -B-j ) (-j 


1 


1 


+mi<i m +< 2 


■)] , 


(77) 
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Table 5. Normalized stress intensity factors for a crack terminating 
at the interface (improved results) 



^+1 


n 2 +i 

y°> 


6 


0 

3.2694 


10 


0 

3.3086 

— = 4 

y l 

2 


1 

3.3149 

3 


1 

3.3420 

= 0.70483 

3 


3 

3.3550 

= 1.29517) 

5 


2 

3.3553 


10 


1 

3.3553 


10 

EXACT 

2 

3.3558 

3.3558 


^+1 


n 2 +i 

r 3 (°) 


6 


0 

13.61 


10 


0 

14.49 

Uo 

25 


0 

15.59 

= 23.077 
y l 

1 


1 

14.89 

2 


1 

16.62 

( 6-j = 0.86935 

2 


2 

17.11 

5 


1 

17.26 

e 2 = 1.13065) 

3 


3 

17.35 


4 


4 

17.35 


5 


10 

17.36 



EXACT 


17.37 


19 






(78) 


k,(b) = 


,+ *l JZ 


g(b) . 


For a uniform crack surface pressure p(r) = p Q and for one material pair 
the plane strain and plane stress results are given in Table 6. The normal- 
ized stress intensity factors shown in the table are defined by 


15,(0) 


k,(0) 

P 0 (b/2) 1 


15, (b) = 


k, (b) 

P (/E7? 


(79a ,b) 


Table 6. Normalized stress intensity factors for a crack ter- 
minating at the interface (normal loading). 



K,(0) 

15, (b) 

Epoxy-Aluminum 

u 2 

— = 23.077 
^1 

Plane strain 
8 1 = 0.6619 

2.7997 

0.8826 

Plane stress . 
8 1 =0.7110 

4.2321 

0.8787 

Aluminum- Epoxy 
Uo 

— = 0.0433 
y l 

Plane strain 
B-, = 0.1752 

0.0981 

1.3421 

Plane stress 
B 1 = 0.1758 

0.0955 

1.3398 


The results shown in Table 6 are obtained by using ten terms in the series 
given in ( 72) and no convergence difficulties were encountered in the calcu- 
lations. 

In conclusion one may note that the technique would be readily appli- 
cable to more complicated problems involving, for example, finite dimensions, 
multiple cuts and general nonsymmetric loading conditions. In this case 
the coupled system of singular integral equations would have Fredholm as well 
as generalized Cauchy kernels and the fundamental functions may again be 
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obtained by using the complex function technique described in this paper. 
Also, one should again emphasize the importance of examining the second (and 
subsequent) roots of the characteristic equations and, if necessary, taking 
them into consideration from the viewpoint of the convergence of the calcula 
tions and the accuracy of the results. 
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